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Abstract —Bridges that are designed to transpose large valleys generally show low stiffness to the lateral bending 
and, therefore, it is not convenient to use the static modeling of the wind actuation. In this context, due to the large 
volume of concrete required for the molding of the pillars, it is fitting to design them with a cross section composed 
of thin walls, in order to confer economy, and with one of the faces open and braced by lintels to guarantee greater 
stiffness to the deformation by bending due to the performance of the shear forces. Such dynamic modeling of the 
pillars activates bending-torsion and presents stiffness to the structure bending divided into two matrices: [J] and 
[S]. Due to the structural stiffness not being bonded in a single matrix, Rayleigh's proportionality in the damping 
matrix assembly [C] does not apply. Thus, the formulation for such an assembly of [C] is proposed in this paper 
via the first three vibration modes of the structure. It models, with examples of application, a bridge with a deck 
backed on three pillars shaped as structural core braced by lintels and of thin-walled sections. Such pillars, with 
axes referenced in the torsional center of each, are inclined off in relation to the global referential of the bridge 
coordinates, in order to apply the decoupling process of the differential equation system that governs the 
phenomenon and the resulting referential transformations. The vibration modes are verified by modeling with the 
ANSYS 2019 R1 software academic version. 

Keywords —Dynamic decoupling, Bridge pillars, Wall panels, Proportional damping, Bending-Torsion. 


I. INTRODUCTION 

In the dynamic analysis of the usual structures, it is 
convenient to model the stiffness of the elements that 
compose them and group them into the unique stiffness 
matrix [K] of the structure in the global coordinate 
reference. As for the pseudo-force due to the inertia that the 
structure exhibits when moving during the dynamic load, it 
is computed in the product of the mass matrix [M] by the 
acceleration vector of the linear and angular displacements 
{d}. Finally, the viscous damping of the structure is 
modeled by a damping matrix [C] proportional to the others, 
where such proportionality is determined by modal analysis 
and the resulting first two vibration modes. Thus, it is 
carried out following the example of Blume [1], and 
Fleming and Romualdi [2]. 

On the other hand, the analysis of the pillar with cross 
section composed of thin walls results in two stiffness 
matrices. The first [7] refers to the portion of the bending 


stiffness resulting from the pillar itself, and the second [S] 
relates to the bending stiffness of the lintels that promote the 
pillar bracing along the height at the open end of the cross 
section. When highlighting the equation of the mentioned 
Wall Panels, Smith and Taranath [3] and also Stamato and 
Mancini [4] are mentioned. In such analyses, the static 
imposition of the loading is done. However, in the dynamic 
analysis, which will be carried out here, the inertial pseudo¬ 
force portion is added to the motion, and so the Partial 
Differential Equations system (PDE's) is constituted by 
three matrices and, to that end, the proportionality procedure 
for the assembly of the damping matrix [C] of the thin wall 
pillars is proposed. For such a proposition, the decoupling 
of the PDE system based on Rosman [5] and necessary 
transformations of referential via Gluck [6] are used. 
Finally, the wind effect is adopted as a function of two 
variables q{x, t ), in terms of the height v of the bridge pillar 
and under variation over time g (t). To that end, the 
positioning of the center of mass is needed for such 
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equation, as inspired by Hart [7], Vaicaitis [8], Laier [9] and 
Awruch [10]. 

The modeling of the wind action in the structures is 
processed, in the majority, from the perspective of the 
imposition of static loads. However, for the case of high- 
section thin-walled pillars (in works of art), lateral rigidity 
is reduced, making dynamic analysis of such a request 
important. This article proposes the systematization of the 
dynamic analysis of such requests in the pillars, thus 
providing a tool for manual calculation or even a mechanism 
of verification of the results obtained by modeling by 
commercial software. 

II. THE WALL PANEL THEORY 

The bridge pillar, considered of cross section made of 
thin and open walls, is conveniently analyzed through the 
Wall Panel Theory. Such technique consists of subdividing 
the core-shaped pillar into several walls and equating the 
balance of forces and moments for the internal forces and 
compatibilizing loads at the intersections of said walls. 
Finally, the differential equations of deflections and 
deformations are used. See Figure 1. 



Fig. 1: Wall panels: (a) Drawing dimension notations, (b) 
Wall referential transformation 

Such equation results in Ordinary Differential 
Equations (ODEs) in terms of the rotation v f , third 
derivative from the displacements v, and the shear forces 
[Vf] generated by external loads. The referred equations are 
written in terms of the center of gravity (CG) referential as: 

-\JlWcG) + [SUV CG } = {V f } ( 1 ) 

where: [7] - Bending stiffness matrix of the Wall Panels; 
[S] - Lintel bending stiffness matrix; { v CG } - Vector of linear 
displacements: v (towards y), co (towards z) and angulars^ 


(rotation around axial axis x)\ {Vf}~ Vector of shear forces 
generated by external loading. 

The matrices [7] and [S], as well as the vector {l^} are found 
in Melo [11]. 


III. DYNAMIC FORMULATION 


The static analysis of the Wall Panels, presented in 
equation (1), is expressed with reference in the center of 
gravity (CG). However, in the dynamic analysis, the inertial 
pseudo-force plot [ M].{i)} is, in terms of mass, distributed 
m at the height H of the pillar, with reference in the center 
of mass (CM). In addition, the mass distribution modeling 
under acceleration z is shown in Figure 2. 



Fig. 2: Distributed mass subject to acceleration z 


The aforementioned pseudo-force is expressed as: 



(2 a) 



(2 b) 


where: m - distributed mass; M - total column mass; H - 
column height. 


For the tridimensional structure, there is: 

0 
m 
0 


[dF}=\M].{ v} = 


m 

0 

0 


0i 

0 

7 


v J 


( 3 ) 


where: I p - Column inertia moment; [M] - Mass matrix in 
the CM. 


Because the dynamic plot expressed in equation 3 is 
referenced in the center of mass (CM) and the static plot in 
the center of gravity (CG), the coordinate transformation of 
the (CM) into the (CG) turns necessary. See Figure 3. 
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Fig. 3: Mass and gravity centers: (a) Positioning , (b) 
Referential transformation of the CM into CG 

The relation of the displacements in the Center of mass 
to the detriment of the Center of gravity is worth: 

V CM — V CG ~ dzcMCG'^CG (4 a ) 

Mcm — g) C g + dycMCG'^cc (4 

0CM — 0CG (4 C ) 

Matrically, the equations (4), expressed as law of 
coordinate transformation of the Center of Mass into the 
Center of Gravity, are expressed as: 

I v cm} = \Jm\'{. v cg) 
v cm\ 

M CM | 

0Cm) 

where: { v CM } - Vector of displacements in the Center of 
mass; {v CG } - Vector of displacemtns in the center of 
gravity; [f M ] - Matrix of linear transformation from CM 
into CG. 

Resulting from the transformation Law expressed in 
equation (5), the Wall Panel movement plot of equation (3) 
is rewritten as: 

{dF CM } = [M].{v cm ] = \T m ] t . [M]. [T m ].{v cg ] (6) 

Deriving equation (1) and adding equation (6), results in 
the ODE"s system that governs the dynamic problem of the 
Wall Panels, expressed by: 

-UUv'cg) + [Si free) + [Ml { Vcg } = {V}} (7) 

where: [M] = [T M ] T . [M]. [f M ], [M] - Matrix of referenced 
mass in the CG, and: V* ext = 2. D ± .x + D 2 \ D 2 = 
-qitfi = -ffr {Vf} = Kxf {71*}; V ext = D ± . x 2 + D 2 . x + 
£ 3 ; D 3 = Q + q 1 .H. 


[J*] by linear transformation by rotation matrix [Re]. Then, 
[J*] is transformed into [J**] which is equivalent to the 
matrix Identity [I] via pre-multiplication and multiplication 
by [J*] -1 / 2 in order to maintain the symmetry of the other 
matrices and guarantee subsequent modal orthogonality and 
be able to use transposed matrices instead of inverse 
matrices. Finally, [M**] and [S**] are diagonalized through 
an iterative process until all terms outside the main 
diagonals of the matrices become of inferior module at a 
determined tolerance s. The set of decoupled PDE's in the 
generalized referential, for the nondamped Wall Panels is 
expressed as: 

-q'j'"(x, t) + Sj. q"(x, t) + nij. qj(x, t) = Vf. (x, t) (8) 

where: Sj is the j' th element of the main diagonal of the 
stiffness matrix of the diagonalized lintel bending [s]; nij is 
the j _th element for the diagonalized mass matrix [m]; and j 
is the j _th vibration mode of the structure. 

In the second procedure, the term of the damping 
[C].{ i>""} is added with the damping matrix [C] constituted 
in terms proportional to the other matrices. Thus, the PDE 
system of viscously damped vibration for the Wall Panels is 
expressed by: 

[MUvcg) + [C].{ v%] + [5].{^' g )- {JUv'cg} 

= (Vr'l ( 9 ) 

where: [ C ] = a M . [M] + a s . [5] — ctj. [/], [C] - Viscous 
damping matrix. 

In the present decoupling, three referential 
transformations are made, where: the first transformation 
consists of diagonalizing the mass matrix [M] referenced in 
the CG by means of diagonalization routine (for example the 
Jacobi Method); in the second transformation, [m*] 
becomes the identity matrix [/] by pre-multiplication and 
multiplication by [ra*] -1 / 2 ; and finally, in the third 
transformation, a looping process is applied for the joint 
diagonalization of [S**] and [/**] up to a certain tolerance 
adopted % and £j. See Flow Chart shown in Figure 4. 

And resulting as a set of decoupled PDE's, expressed as: 

qj(x, t) + Cj. q""(x, t) + Sj. q"(x, t) -jj. q""(x } t) 

= V f °j(x, t) (10) 


1 ^ ^ Z CMCG 

0 1 dycMCG 

0 0 1 



(5) 


IV. DYNAMIC DECOUPLING 

The first procedure is to rewrite the system of ODEs 
(equation 7) in the main frame of the bending pillar stiffness 
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Fig. 4: Dynamic decoupling flowchart via Procedure 2 


V. MODAL ANALYSIS 

The ODE is written in the space v and temporal in t by 
means of the PDE's resolution by the Method of Modal 
Superposition for the nondamped vibration as: 

jj.u""(x) — Sj.u"(x) — cof.u(x ) = 0 (11a) 

g(t) — o)j.g(t) — 0 (116) 

where: q(x,t) = u(x).g(tj, C 0 j is the j _th vibration 
frequency of the structure. 

The analysis of equation 11 (a) is carried out in order to 
determine the vibration modes and correlate Sj and jj. Then, 
u(x) = e a - x is applied and the characteristic equation is 
expressed by: 

jj. a 4 - Sj. a 2 - cof = 0 (12) 


Where the roots of the characteristic equation, equation 
(12), aq and a 2 , and the relation jj and Sj in terms of the j _th 
vibration frequency of the structure expressed by: 



(13 a) 



(13 b) 


Jj 


Sj.a 2 + cof 


(13 c) 


The solution of the displacements u(x) is expressed as: 

u(x) = C^sinh^a^x) + C 2 -cosh(a 1 .x ) + C 3 .sin(a 2 .x ) 
+ C 4 .cos(a 2 .x ) (14) 


Adopting the nondimensionalization of the roots a into 
A, via A 1 = a 1 .H and A 2 = a 2 .H , roots A 1 and A 2 are 
correlated by: 

A 2 -A 2 = A 2 (15 a) 



(15 b) 


Also, by applying the nondimensionalized roots and 
a 2 in the form of A 1 and A 2 and multiplying the latter, the 
frequencies of the structure are expressed as: 


' H 

Ai- A 2 


( °j =( °j-- H 2 


CO; 


(16 a) 
(16 b) 


The pillar analyzed here is of the fixed-in-the-base type, 
therefore, with displacement u (x = 0) = 0 null and rotation 
u'(x = 0) = 0 also null. At the top, the end is free and, at a 
subsequent moment, bending and stress are shear and null, 
thus applying: u"(x = H) = u’"(x = H) = 0, 

respectively. Applying such boundary conditions in 
equation (14), after due nondimensionalization, the linear 
equation system results in: 

1 0 1 0 - 

0 2-1 0 A 2 

Al .coshOO Al-sinhOO -0% -cos(A 2 ) -k 2 2 .sin0 2 ) 
^l.sinhOO Al-coshUO .sin(A 2 ) -\\.cosOi)- 
= { 0 } 

(17) 


The nontrivial solution of the system, shown in equation 
(17), leads to the transcendental equation: 

f.cosh^Afj.cos^Af) + f 2 .sinh(Af).sin(Af) = —1 (18) 
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with: f t = 


2.A\.A\ 


a x .a\ + a\.a 2 


h = 


A\. A\ - A\. A\ 
X\. A 2 + A}. A 2 


Highlighting that, the period will be: 



2.7T _ ^ H 2 


CO 


'J 
2. 7T 

A i- A 2 




(19 a) 
(19 b) 


In Figure 5, the change in the elastic line in the first four 
vibration modes is presented for A 2 = 0. Figure 6 shows the 
ANSYS modeling, evidencing some vibration modes of the 
column with open cross section and thin walls. 



— 1st vibration mode 

— 2nd vibration mode 

— 3rd vibration mode 
4th vibration mode 


Fig. 5: The first four vibration modes of the Pillar m with 
fixed-at-the-base core 

Now, through the relation between the 
nondimensionalized roots A 1 and A 2 , constant in equation 15 
(a), and in the substitution of the harmonic and hyperbolic 
terms for the Taylor Series with 11 terms, the vibration 
modes j = 1 to 7 are expressed by A = 0, see Table 1. 


Table. 1: Vibration modes of the Core Column by A 2 = 0 



Mode 

Mode 

Mode 

Mode 

Mode 

Mode 

Mode 


i = 1 

i = 2 

i = 3 

i- 4 

i- 5 

i- 6 

i = 7 

A\ 

1.87510 

4.69409 

7.85483 

9.08911 

10.02092 

11.68901 

14.62618 

A'l 

1.87510 

4.69409 

7.85483 

9.08911 

10.02092 

11.68901 

14.62618 

co\ 

3.51600 22.03448 61.69835 84.61192105.97714136.63295213.92514 

T\ 

1.78703 

0.28515 

0.10184 

0.07606 

0.05929 

0.04599 

0.02937 


A: Modal 

Total Deformation 2 
Type: Total Deformation 
Frequency: 0,28489 Hz 
Unit: m 

11/04/2019 00:01 

0.0022751 Max 

0,0020223 
0,0017695 
-| 0,0015168 
Pi 0,001264 
0,0010112 
—I 0,00075838 
—| 0,00050559 
0,00025279 
OMin 


m 


ACADEMIC 


Y 

X 


A: Modal 

Total Deformation 5 
Type: Total Deformation 
Frequency: 1,7551 Hz 
Unit: m 

11/04/2019 00:03 

m 0.0022886 Max 

0,0020343 


— 0,0015257 

— 0,0012714 

— 0,0010171 
—| 0,00076285 

0,00050857 
0,00025428 
OMin 


A: Modal 

Total Deformation 14 
Type: Total Deformation 
Frequency: 7,0834 Hz 
Unit: m 

11/04/2019 00:03 

m 0,0032475 Max 

0,0028866 
0,0025258 
—I 0,002165 

— 0,0018042 

— 0,0014433 
—| 0,0010825 

0,00072166 
0,00036083 
OMin 


Fig. 6: Vibration modes by modeling with the ANSYS 2019 
R1 software academic version 


VI. PROPORTIONAL DAMPING 

Procedure 2 of diagonalization of the PDE's system 
presented in equation(9) results in: 

Cj = a M . 1 + a s .Sj — aj.jj ( 20 ) 

The analysis of the PDE shown in equation (10) and 
subsequent application of u(x) = e ax and g(t) = e a)t , 
results as a characteristic equation of the damped free 
vibration: 

a) 2 + ( Cj.a 4 ).co + ( Sj.a 2 — jfa 4 ) = 0 ( 21 ) 

The resolution of the characteristic equation, equation 
(21), with the imposition of critical damping and subsequent 
u of the damping ratio <(y implies cf equals to: 

cf = 2. <fj. Mj (22) 

where: f j - damping ratio of the J th vibration mode. 

Equations 13 (c) and 15 (b) with due 

nondimensionalization of the roots A 1 = a 1 .H and A 2 = 
a 2 .H , result in Sj and jj in terms of the frequency Ci)j, as: 
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007 


j] (A\-A\.A 2 ) (23a) 

Sj=A 2 .jj (23 b) 

Finally, since it is only desired to establish proportions 
of the matrices [M], [/] and [S] that form the damping 
matrix [C], the application of equations (22) and (23) in the 
condition expressed in equation (20) and using the 
frequency oo] instead of ooj , since it is only desired to 
establish proportions of the matrices results in : 


_ &m a s- s j a j-jj 

' j ~ Zoo] + 2 Too] Zoo] 


(24) 


Imposing the first three vibration modes j = 1,2 and 3 
in equation (24), a system of equations is built, whose 
solution is: 


% : 


a s ■■ 


. 2.K-<fi0 2 -s 3 ~js-s 2 ) + (Oz-flVs-Si -ji.Ss) + <o$.Z 3 (ji.s 2 -; 2 -Si)] 


j 1 - (s 2 - S 3 ) + j 2 . (s 3 - Si) + j 3 . (s x - s 2 ) 


. 2 - hA-A(A -) 2 ) + a*2-&(/i - A) + ^3-^3 (A - A)] 

A- fe - S3) +72- (S3 - Si) +A- (5! - S 2 ) 


(25 a) 
(25 b) 


= 2 - hAjife - S 2 ) + (A-AGi - S 3 ) + qjgjsfe - %)] f25r"> 

A- (s 2 - s 3 ) + A- (s 3 - sj + A- (si - s 2 ) ^ ^ 

In percentual terms there are: 


Mm - 

Ms = 

My = 


a M 


a M + % + a y 

a s 


a M + % + a 


7 


a 


y 


a M + a 5 + a y 


(26 a) 
(26 b) 
(26 c) 


a) Center of Mass Positioning 

For the reinforced concrete, cross, and thin-walled open 
section, the determination of the center of mass coordinates 
(CM) is made necessary. To this end, the following 
definitions are used: 

Yi=i x CMf m i 


Xtm — 


SI Umt 

I.i=i y cm r m i 


y cm 


i 27 (J) 
(27 b) 


Z?=i m t 

where: x* CM . andy£ M - coordinates of the canonic center of 
mass i in relation to a generic referential x* and y* defined 
in the cross section; m t - canonic element mass i\n - 
number of canonic masses that formed in core; andxJ M and 
y] M - the center of mass coordinates of the cross section in 
core shape in relation to referential x* and y*. 


The adoption of the subdivision presents in Figure 7, the 
concrete specific mass p cs = 2400 kg/m steelp^- = 
7860 kg/m 3 results in: 

X CM^ m l — x CM 2 - m 2 ~ X CM 3 - m 3 ~ X CM 4 -^4 + X CM 5 -M 5 


X CM ~ ' 


m 1 — m 2 — m 3 — M 4 + M s 


(28 a) 


* Vcm^i ~ y*cM 2 -^2 - y*cM 3 - m s - y*cM 4 -M 4 + y* CMs -M s 

- mi - m2 - TO3 - M4 + Ms - (28fc) 


B V by 


1 : 


X CM ^ CM 

r 

t y'cM 


Bl a 'y 


t ycM 

i®E |l a2v 

.• ! ♦ t. — •> 

b z 1 1 '2 x * 


4 XCMl - CMX 

y'cM, 


(-A 




Xct *2 CM 2 

« 1 
A 

ycM 2 


(-) 




•I *• 


0 


b z + 


L I 

^ jL J y'cM, 


r \ 


F.“ 0 -,-o- —»Tool j 7, 

o 


(-) b v h, 


y y cp 
3 


!::! x cn< erf 

Bl t 


ycM 4 


( + ) By by 


-1 £o_o 0000 gtt <L»o_oooj X t, 

1. 1 h 


ST '2 
a iy 


r.7»|y 

l CH " f 


u-i—J 

B z 


t § 1 U2y 

. . - 


Fig. 7: Geometry of decomposition into canonic format 
masses 

Finally, the distances between the centers of mass (CM) 
and that of gravity (CG) are valid, according to what is 
shown in Figure8: 

y* 



•••••••••••••••••• 

:! c 

Ci 

G < 

M 

^ dycMCG 



11 








“ZCMCC “ X CM X CG 

Fig. 8: Positioning of the CG and CM in the cross section 
d y C MCG = y*0M - y*CG (29 a) 


*’■ZCMCG X CM X CG 


(29 b) 
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VII. REFERENTIAL TRANSFORMATION 

a) Transfer from the center of gravity CG to the torsional 
center D 

According to the compatibilization of the rotations (p D 
(in the torsional center) and f CG (in the center of gravity), 
the displacements o) D via co CG , and projection of the distance 
dy CG D from the center of gravity to the torsional center 
parallel to z axis, procedures are the same forv D with v CG 
and d ZCGD . See Figure9. 



Fig. 9: Transformation of the referential in the CG for D: 
(a) translation, (b) rotation 


The matricial format of the displacement balance is: 


{Vcg } — I Td\- { v d) 


' v cg' 


1 0 d ZcGD 


-Vd' 

MCG 
$CG; 

— 

0 1 dycGD 

0 0 1 


(l) D 


(30 a) 
(30 b) 


b) Transfer from the torsional center D for the source of the 
coordinates 0 5 

Analogously as in (a) of this item, and considering the 
rotation (3 between the referential ( yDL> z DL> x DL ) in the 
torsional center and the global referential ( y g ,z g ,x g ), see 
Figure 10, there is: 



Fig. 10: Transformation of the referential in the torsional 
center D for the source 0 S . 

v D .cosjl + co D .sin(l — v g — d ZDOs .4> g (31 a) 

—v D .sin(l + cjl> d .cos(J = co g + d yDOs . f g (31 b ) 


<t>D = 4>g 

Matrically, the balance is expressed as: 
{v D } = [T s ]~\[T G l{v g } 

10- 
0 1 



cos(3 sin(3 0 
—sinfl cosfl 0 
0 0 1J 


L0 0 


(31c) 


(32 a) 



(32 b ) 

Combining equations (30) and (32), direct transfer from 


CG to O s is as follows: 


{v g } = [T s ]- 1 .[T G ].[T D ]-\{v cc } 


(33) 


VIII. RESULTS AND DISCUSSIONS 


In this section, a bridge with a deck backed over three 
pillars, each with a cross section composed of thin walls will 
be analyzed. See Figure 11.All in compliance with the 
Technical Standards NBR 7188 [12] and DIN 1055 [13]. 


The material used is reinforced concrete of resistance 
class C - 40, complying with NBR 6118 [14]. Therefore, the 
modules are of longitudinal elasticity E = 
3.54175 x 10 7 kN/m 2 and of cross elasticity G = 
1.47573 x 10 7 kN/m 2 , as well as Poisson's coefficient v = 
0.20. On the basis of the loading conditions of the bridge, 
two scenarios are analyzed: in the first scenario, the sole 
actuation of the wind on the structure (See Example 1) and, 
in the second scenario, the traffic performance with a focus 
on braking (See Example 2). 


. In this first analysis, the bridge is subject to 
the sole actuation of the lateral wind in the form of three 
loadings along the height of the pillar, namely: q 1 = 
20 kN/m, q 2 = 30 kN/m and Q = 10 kN. See Figure 12. 




Fig. 11: Bridge backed over three thin-walled pillars: (a) 
lateral view, (b) drawing configuration of the pillars 
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Fig. 12: Angles formed between the external shear 
forcesV extp ^ V extp ^, V ext p 3 and axes y DL defined in the 

torsional centers of each of the three columns 

Vexty ~ — 2~f ]' % ~ Qly ^0 

Qly = Ql- cos (@ d ) > 4ly — (fo- cos (^ d ) f Qy 
= Q.cos(O d ) 

q lz = <7i-sin(0 D ); q 2z = q 2 . sin(0 D ) ; Q z = Q. sin(0 D ) 
V'ext, = ~Y ^- x2 ~ + C<?Z + q lz .H) 


The shear force, over time, will be governed by a 
harmonic sinusoidal function V g (t) = sin(co.t), with 
frequency co = 3 rad/s. In addition, the PDE system of 
equation (9) is expressed for the pillars and referenced in the 
center of gravity as: 
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The same PDE system, with the proportionality 
coefficients of the matrices [M], [7] and [S] to compose the 
viscous damping matrix [C], is expressed as: a M = 
0; a s = -3.80 %10 16 ; - -3.80 x 10 16 ; = 

0; q. s = 0.5; \x ] — 0.5. 


And the set of decoupled PDE's presented in equation 
(10) are: 

-q^"(x,i) + 0.03385. q^x, t) + 0.49519. q""(x, t) 

- 0.00962 .q"(x,t) = Vf x (x, t) 

-q 2 "(x t t) + 0.04038. q 2 (x, t) - 0.13363 .q' 2 '(x t t) 

- 1.26726 .q 2 (x,t) = V? 2 (x,t) 

-q'{"(x, t) + 0.00003 .q 3 (x, t) + 0.50243 .q 2 "(x, t) 

+ 0.00498. q'f(x,f) = V^(x, t) 


(34) 


In the set of EDP’s presented in equations (34), the first 
mode of vibration (mobilizing the lintels) is verified by 
crossing the equation (15 b) and the Figure 13. 
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A: Pilar em C com Linteis - C-40 

Total Deformation 4 
Type: Total Deformation 
Frequency: 1,7528 Hz 
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A: Pilar em C com Linteis - C-40 

Total Deformation 6 
Type: Total Deformation 
Frequency: 4,4789 Hz 
Unit: m 
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Fig. 13.‘Column C in concrete C - 40, vibration modes 
activating the lintels for flexion: (a) 1st mode (co 1 = 
1,91493 rad/s), (b) 2nd mode and (c) 3rd mode via 
modeling in ANSYS 

so, with the parameters Sj = —0,00962 and jj = — 1, and 
through equation (15 b) the parameter is determined A = 

H-JV/ = 100 [m], 0,0096 y_ 1 = 9,81, using Table 

2 with the value immediately above X, then: of* = 
2,54634 rad/s (for A = 10). 

Table. 2: Vibration modes of the Core Column by A = 10 



Mode 

i = 1 

Mode 

i = 2 

Mode 

i = 3 

Mode 

i = 4 

Mode 

i = 5 

Mode 

i = 6 

Mode 

i = 7 

A± 

3,25746 

5,50391 

8,42839 

9,64266 

10,43138 

11,87014 

14,61565 

^2 

0,78169 

4,50478 

7,81267 

9,10938 

9,94051 

11,44116 

14,26945 

<*>1 

2,54634 

24,79389 

65,84820 

87,83870 

103,69321 

135,80821 

208,55730 

rp * 

1 i 

2,46754 

0,25342 

0,09542 

0,07153 

0,06059 

0,04627 

0,03013 
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Finally, in Figure 13 (a) the first vibration mode of a 
single C-shaped pillar and braced with lintels is presented 
under the same distribution as this example, hence the 
reason for comparison, with an approximation of 2.24 
%.The frequency adjustment determined via TMC is used 
by the coefficient a = Vl 3’ whose functionality is based 
on Diziewolski [15]. 
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The shear force Vf functions in the initial referential are 
presented in Table 3, as well as Vf in the generalized 
referential, and both for the columns P l5 P 2 and P 3 . 


Table. 3: Functions of the shear force per Column 



Column^! 

ColumnP 2 

ColumnP 3 


89.77 .x 2 

-48.87 .x 2 

-147.66.x 2 

V fl 

+ 11969.20 .x 

-6515.63.x 

- 19687.64.x 


- 1196926.12 

+ 651565.87 

+ 1968773.75 


120.17.x 2 

-141.82 .x 2 

26.41 .x 2 

% 

+ 16023.05.x 

- 18.908.90.x 

+ 3520.92 .x 


- 1602313.28 

+ 1809899.76 

- 352093.65 


332.19.x 2 

-180.83 .x 2 

-546.40.x 2 


+ 44291.71.x 

-24110.90.x 

- 72853.59.x 


-4429193.40 

+ 2411102.23 

+ 7285395.07 


0.18.x 2 

-0.10.x 2 

-0.30.x 2 

V fi 

+ 24.43.x 

- 13.30.x 

-40.17.x 


- 2442.55 

+ 1329.69 

+ 401750.64 


0.27.x 2 

-0.32.x 2 

0.06.x 2 

v° 

V f2 

+ 35.85 .x 

-42.31.x 

+ 7.88 .x 


- 3585.47 

+ 4231.26 

- 787.97 


0.59.x 2 

-0.32.x 2 

-0.96.x 2 

17° 

V h 

+ 78.21.x 

-42.57.x 

- 128.64.x 


-7820.51 

+ 4257.22 

+ 12863.63 


For the first vibration mode, the displacement function, 
u(pc) and g(t ), and the generalized and initial referential 


are, respectively: 

u ± (x) = 0.36275 .{cosh( 0.01983 .x) 

- cos(0.01713 .x) 

+ 0.92535. [- sinh(0.01983 .x) 
+ 2.41421.sm(0.01713 .x)]} 


0i(t) = 4.543.10 10 . sin(3.476.10“ 3 . t) 
- 8.928.10 8 .sm(3.t) 


u^x) = 6.098.10 -4 . {cos/i(1.847.10 -3 .x) 

— cos(7.65.10 -4 .x) 

+ 1.556.10" 3 . [- sinh(1.847.10~ 3 . x) 
+ 4.048.10" 3 . sin(7. 65.10" 4 . x)]} 

flfi(t) = 7.637.10 7 . sin(5.843.10" 6 . t) 

- 1.501.10 6 .sm(3.t) 

Then, returning to the initial referential, the vector of 
displacements { v CG } for the mode of larger deformations is 
written as: 

{Veal 1 = {<?} = 

1 

where: [® R ] = [Re]. [/*]”. [<p ± ]. [<p 2 \. [<p 3 ]. (...). [<p n -il [<Pnl 

And, finally, the rotation function <p CG (x) is expressed 
as: 


0 CG (x) = u^xfg^t) + u 2 (x).g 2 (t) + u 3 (xfg 3 (t ) 


The relation between rotation 4>cgM an d the 
Bimomentf?(x), originating from the Bending-Torsion 
Theory yields: 


B(x) = E.I 0J . 0" G (*) = E.I 0) . 


rf 2 [^ccW] 

dx 2 


After applying the referential transformations, we have 
the Bimoment in the Base and in the Top for the CG, D and 
O s . See Table 4, as well as the rotation in Table 5. 


Table. 4: Bimoment values per Pillar in the CG and D, and 
the set of pillars in O s 


(x 10 7 ) 

kN.m 2 


B(x = 0) 
B(x 

= 50 m) 

B O = H) 


In the Center of In the Torsional 

Gravity Center q 

CG\ CG 2 GG 3 D-i D 2 TJ 3 

-3.99422 2.17432 6.56993 -3.99422 2.17432 6.56993 

Idem 

-4.00790 2.18176 6.59242 -4.00790 2.18176 6.59242 per 

Pillar. 

-4.04993 2.20464 6.66156 -4.04993 2.20464 6.66156 


Table. 5: Rotation values per Pillar in 
pillars in O s 

the CG and set of 

In the Center of Gravity 

rnd 

o s 

CG t 

cg 2 

cg 3 

II 

o 

o 

0 

0 

Idem 

0 (x = 50m) - 57.42669 

31.26114 

94.45877 

per 

00 = H) - 233.82458 

127.28615 

384.60827 

Pillar. 
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EXAMPLE 2: In this second analysis, there is the braking 
force of the traffic and lateral cargo at the top. Thus, the 
shear force under the configuration shown in Figure 14 is 
presented as: 



Fig. 14: Angles formed between the external shear forces 
activated by braking in each Pillar and the axes y DL in each 
of the three columns 


Adopting the same frequency co = 3 rad/s for the shear 
force, the values of the Top and Base Bimoments are 
expressed in Table 6. 

Tab. 6: Values of Bimoment for the actuation of the 
braking force of the vehicles 


O 10 5 ) 

kN. m 2 - 


In the Center of 
Gravity 


In the Torsional 
Center 


O x 


CG'i CGo 


CGo 


D i 


D 7 


Do 


50 = 0) -4.96373 2.70209 8.16464 -4.96373 2.70209 8.16464 


-4.98075 2.71135 8.19264 -4.98075 2.71135 8.19264 


B O 
= 50m) 

B(x = H) - 5.03302 2.73980 8.27860 - 5.03302 2.73980 8.27860 


Idem 

per 

Pillar. 


where the loadings are: q x — q 2 = 0 kN /m, Q = 10 kN ; 
the Top Bimoment B H = 10 kN.m 2 ; and the torsional 
moment at the top is null M tH = 0 kN. m. 


IX. CONCLUSION 

In this paper the implementation of the dynamic 
plots of acceleration and damping in the Wall Panel Theory 
was carried out, in order to promote the dynamic analysis of 
bridge pillars under the effect of Bending-Torsion and to 
provide the noncoincidence of the center of gravity with the 
torsional center, as a result of the geometry in thin-walled 
cross section. Moreover, the PDE system that governs the 
vibration of said Wall Panels has two bending stiffness 
matrices, being: [7] for the pillar itself and [ S ] for the lintels 
that promote bracing of said pillar. Thus, the Rayleigh 
proportional [C] damping matrix assembly considers only a 
single stiffness matrix. At this point, an unpublished 
formulation is made for such a composition. Two 
procedures of diagonalization of the PDE system and 
transformations between the reference points of the center 
of gravity (CG), the center of torsion ( D ) and a global origin 
(O s ) are discussed. 


Such analyses encourage the of high columns in 
reinforced concrete subject of wind actuation, due to the 
large volume of concrete required for the molding of the 
pillars. It is fitting to design them with a cross section 
composed of thin walls, in order to confer economy, and 
with one of the faces open and braced by lintels. 
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